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CN ■ ABSTRACT 

>>■ 

f^ , Context. That the rotation curves of spiral galaxies are generally flat in the outer regions is commonly considered an indication that 

galaxy disks are embedded in quasi-isothermal halos. In practice, disk-halo decompositions of galaxy rotation curves are performed 

in a parametric way by modeling the halo force contribution by means of expressions describing approximately the properties of the 

regular isothermal sphere or other spherical density distributions suggested by cosmological simulations. 

Aims. We construct self-consistent models of nonspherical isothermal halos embedding a zero-thickness disk, by assuming that the 

halo distribution function is a Maxwellian. The method developed here can be used to study other physically-based choices for the 

halo distribution function and the case of a disk accompanied by a bulge. 

Methods. The construction is performed by means of an iterative procedure that generalizes a method introduced in the past. In a 

preliminary investigation we note the existence of a fine tuning between the scalelengths R^ and h, respectively characterizing the 

\^ , rise of the rotation curve and the luminosity profile of the disk, which surprisingly applies to both high surface brightness and low 

surface brightness galaxies. This empirical correlation identifies a much stronger conspiracy than the one required by the smoothness 

and flatness of the rotation curve {disk-halo conspiracy). 

Results. The self-consistent models are found to be characterized by smooth and flat rotation curves for very different disk-to-halo 

Q ■ mass ratios, hence suggesting that conspiracy is not as dramatic as often imagined. For a typical rotation curve, with asymptotically 

$H ' flat rotation curve at V™ (the precise value of which can also be treated as a free parameter), and a typical density profile of the 

ri ' disk, self-consistent models are characterized by two dimensionless parameters, which correspond to the dimensional scales (the disk 

j^ I mass-to-light ratio M/L and the halo central density) of standard disk-halo decompositions. We show that if the rotation curve is 

decomposed by means of our self-consistent models, the disk-halo degeneracy is removed and typical rotation curves are fitted by 

models that are below the maximum-disk prescription. Similar results are obtained from a study of NGC 3198. Finally, we quantify 

the flattening of the spheroidal halo, which is significant, especially on the scale of the visible disk ("abridged"). 

> ■ 
^^ ' Key words, galaxies: spiral - galaxies: structure - galaxies: halos - galaxies: kinematics and dynamics - galaxies: individual NGC 

m ' 3198 

rn 

(y-^ 1. Introduction photometric profile and the corresponding velocity data-points 

f~*) can be modeled in many different ways, either by heavy or by 

(^ The study of the rotation curve of spiral galaxies (from the case light disks, through an ill-determined mass-to-light ratio M/L 

,-H , of M31, van de Hulst et al. 1957, to the systematic investigation for the disk. In general, the available data are only able to set 

•• ■ of nearby galaxies in the 70's) has led to convincing evidence for an upper limit to such M/L, identifying the so-called maximum- 

.!^ ' the presence of dark matter halos (e.g., see Roberts 1976). To disk solution, for the model that maximizes the contribution of 

^ this discovery, HI rotation curves, which often sample regions visible matter Typically, the maximum-disk solution is able to fit 

;_, well beyond the bright optical disk, played a decisive role (van the rotation curve out to the optical radius (Kalnajs 1983; Kent 

C3 , Albada et al. 1985; van Albada & Sancisi 1986). 1986). 

In this context, it has been suggested that the existence of a It was immediately realized that the flatness of the outer ro- 

halo of dark matter introduces a problem of fine tuning or "con- tation curves suggests that the dark halo embedding the visible 

spiracy" (van Albada & Sancisi 1986). The coexisting dark and disk might have a structure similar to that of the classical isother- 

visible components, while expected to have different formation mal sphere, because outside the central core such sphere is char- 

and evolution histories (even if naturally coupled by their mutual acterized by a density distribution declining approximately as 

gravitational interaction), must always cooperate to produce a p ~ l/r^. Curiously, the ubiquitous presence of gravitational 

gravitational field inversely proportional to the radius, so that the fields consistent with a flat rotation curve also came from studies 

resulting rotation curve remains flat from the inner regions (gen- of elliptical galaxies, not only in the nearby universe (for NGC 

erally thought to be dominated by the visible matter) out to large 4472, see lower left frame of Fig. 7 of Saglia et al. 1992 and 

radii (where the field is dominated by the diffuse dark halo). This Bertin 1993; for IC 2006, see Franx et al. 1994; for a more recent 

puzzling feature is accompanied by a serious problem of degen- systematic study, see Gerhard et al. 2001), but also for ellipticals 

eracy. In fact, from a detailed study of the so-called disk-halo at cosmologically significant distances (e.g., see Barnabe et al. 

decomposition, it has been shown that, in general, the available 2009, and references therein). 
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Of course, deviations from the paradigmatic case of galaxies 
such as NGC 3 198 are not uncommon (e.g., see Casertano & van 
Gorkom 1991). But it has been pointed out that such deviations 
are generally associated with well-defined features present in the 
distribution of visible matter or the presence of non-circular mo- 
tions (Sancisi 2004). On the other hand, pure disks characterized 
by the absence of specific photometric features are quite com- 
mon (see also Kautsch 2009). 

The literature on the general topic of the distribution of dark 
matter in spiral galaxies is enormous and deals with a variety 
of important aspects of the problem. In recent years these stud- 
ies have generally focused on the requirements of current cos- 
mological scenarios. In this Introduction it would be impossible 
to do justice to all the relevant studies in this general research 
area. However, it appears that one simple aspect of the problem 
has gone practically unnoticed (but see some discussion given 
by van Albada & Sancisi 1986). This is the issue of the mutual 
interactions between halo and visible matter in determining the 
self-consistent equilibrium of the basic state. In their review, van 
Albada & Sancisi (1986) addressed this point, with the help of 
N-body simulations, but did not carry out a systematic investi- 
gation. 

Here we wish to reconsider this problem, starting from an ex- 
tremely simple physical (and somewhat academic) description of 
the dark halo. We will thus address the problem of determining 
the nonspherical structure of a dark halo, assumed to be perfectly 
isothermal, in the presence of a given disk. We will show that 
even if very simple, and thus lacking some important features 
present in real objects (for example, the presence and role of gas 
are not considered explicitly in this paper, although they could 
be easily incorporated), this study leads to some interesting and 
unexpected results. 

The search for equilibrium solutions of the Vlasov-Poisson 
equations, where the mean potential is generated by the distri- 
bution of stars populating the system, with the possible contri- 
bution of an "external" density distribution (e.g., associated with 
a co-existing dark matter halo), is a classical problem of stellar 
dynamics, with important analogues in the study of collision- 
less plasmas. In this paper we follow the approach based on the 
Jeans theorem to deal with the halo matter, while the disk of 
visible matter is taken to provide an external field. In this ap- 
proach, the starting point is the identification of an appropriate 
form for the one-particle distribution function in phase space; 
in our case this is the isothermal assumption for the dark halo. 
The adopted distribution function then generates a O-dependent 
density in physical space, thus leading, inserted as a source in 
the Poisson equation, to a nonlinear problem in €). Here O rep- 
resents the gravitational potential. In the past, this method has 
been used frequently and has led to important insights into the 
dynamics of galaxies and other stellar systems. Yet the study of 
self-consistent equilibria in the presence of external fields that 
break the natural symmetry associated with the one-component 
problem is only rarely considered (see Bertin & Varri 2008). 
Therefore, besides the natural application to the discussion of 
the structure of dark halos around galaxies, this paper has some 
interest of its own, as a study in stellar dynamics. 

The structure of the paper is the following. In Sect. 2 we 
comment on the so-called disk-halo conspiracy. Here we iden- 
tify an empirical correlation which brings out a new and stronger 
conspiracy than the one associated with the flatness of rota- 
tion curves. In Sect. 3 we formulate the self-consistent prob- 
lem of constructing an isothermal spheroid in the presence of 
a zero-thickness disk. The actual procedure to calculate the self- 
consistent models is reported in Appendix A. The properties of 



self-consistent isothermal halos are then described in Sect. 4. 
Then in Sect. 5, by modeling first a fiducial rotation curve and 
later by modeling the specific case of NGC 3198, we show that 
self-consistency actually removes the degeneracy that charac- 
terizes standard parametric disk-halo decompositions (with ad- 
ditional comments given separately in Appendix B). Finally, 
Sect. 6 summarizes the main results of the paper. This article 
is all based on the mathematically simple case of zero-thickness 
disks. A follow-up paper will address the corresponding study 
of disks with finite thickness. 



2. A comment on the disk-halo "conspiracy" 

The rotation curve V{R) of many spiral galaxies is featureless, 
regular, and smooth: after a central region of approximately lin- 
ear growth, the rotation velocity reaches a constant value Voo 
and remains flat out to large radii, often well beyond the bright 
optical disk. [But, of course, this is meant to be only a zero-th or- 
der description, because even flat rotaton curves are not entirely 
featureless (Sancisi & van Albada 1987).] Since the underlying 
gravitational field results from the sum of the contributions of 
different mass components of the system (the stellar disk, the 
bulge, the gaseous disk, and the dark matter halo), which gen- 
erally dominate in different regions, the smoothness and lack of 
features of the rotation curve is unexpected and poses a physical 
problem, which is commonly referred to as the disk-halo con- 
spiracy. 

When such conspiracy was noted (see van Albada & Sancisi 
1986) most of the emphasis was placed on the flatness of the 
rotation curves. In reality, it has long been noted that standard 
parametric disk-halo decompositions, based on the superposi- 
tion of the fields produced by the visible matter and by a spheri- 
cal halo, are subject to a high level of degeneracy. Equally viable 
fits to the rotation curves are obtained with very different models 
(i.e. very different values of the mass-to-light ratio of the stellar 
disk, assumed to be constant within the galaxy), ranging from 
some characterized by a very light disk and a dominant halo all 
the way in, to others, close to the so-called maximum disk so- 
lution, where the support to the rotation curve inside the optical 
disk is basically due to the visible matter alone, while the halo 
contribution becomes important and eventually dominates only 
outside. In a sense, the existence of such degeneracy actually 
undermines, or at least softens, the problem of conspiracy ap- 
parently raised by the flatness of rotation curves. 

Of course, there are several constraints that can be used to re- 
duce the level of degeneracy present. In addition to those posed 
by stellar population analyses, we may recall the constraints 
posed by dynamical stability in the context of spiral structure 
(Ostriker & Peebles 1973; Bertin & Lin 1996; see also Widrow 
et al. 2008; and references therein) or warp dynamics (e.g., see 
Bertin & Mark 1980). 

To give a concrete example of the procedure adopted in 
standard parametric disk-halo decompositions, we may recall 
the case of a purely exponential disk with the adoption of the 
model of a classical pseudo-isothermal sphere for the halo con- 
tribution. In this case the rotation curve is decomposed into 

^mod ^ DM ^ ^ D- "^^^ 



vliO = vi^f 



'«(f)-«(l)-"(l)Mf) 



(1) 



where ^ - R/his the galactocentric cylindrical radius in units of 
the exponential scalelength and /? is a dimensionless parameter 
(see Eq. ( fT6b ) which measures the weight of the stellar disk. In 
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addition, /q, Kq, h, and Ki denote the standard modified Bessel 
functions. The halo contribution is described by 



vIm(^) = vi 



1 



iiarctan(^^ 



(2) 



where a is a dimensionless measure of the central dark matter 
density (see Eq. (fT9]l). In this example, the degeneracy occurs in 
the {a,/3) parameter space (see Sect.jSj. 

In the following subsection we will show the existence of a 
much stronger reason for invoking a problem of conspiracy, as- 
sociated with the fine tuning between the length scale that char- 
acterizes the inner growth of the rotation curve and the natural 
scale length of the visible disk (in the absence of significant con- 
tributions from a bulge or from gas, this would be the exponen- 
tial length). This tuning is a direct indication of a subtle interplay 
between the density distribution of dark and visible matter If 
the two length scales were not sufficiently close to the observed 
correlation, the degeneracy of the standard parametric disk-halo 
decompositions would not take place. 

2.1. An important empirical correlation 

This subsection will provide a useful empirical background to 
the physically simple dynamical models presented in this paper. 
Here we start by identifying one important aspect of the prob- 
lem of the uniformity of the observed characteristics of rotation 
curves. 

For simplicity, we refer to a sample of galaxies for which the 
contribution of the bulge should be negligible. We first consider 
disk galaxies belonging to the Ursa Major cluster, for which 
good-quality photometric profiles (Tully et al. (1996)) and rota- 
tion curves (Verheijen & Sancisi (2001)) are available. Without 
actually aiming at producing a complete sample in any statisti- 
cally significant sense, which would go well beyond the goals 
of this paper, we then add a small number of other galaxies, 
with data taken from Kent (1987), de Blok & McGaugh (1997) 
(with photometries from de Blok et al. (1995)), and Swaters et al. 
(2000) (with photometries from de Blok et al. (1995) and from 
McGaugh & Bothun (1994)), so as to increase the size of the 
sample. The only criterion used to select the galaxies from the 
source papers is the shape of their photometric profiles. Galaxies 
with clear deviations from an exponential profile in the central 
parts, often related to the presence of a bulge, are discarded. 
Therefore, for these approximately exponential stellar disks, any 
difference in the shape of the rotation curves should be inter- 
preted only in terms of the different relative weights of the disk 
and of the dark component (under the simplifying assumption 
that the mass-to-light ratio is constant within each galaxy disk). 
The final sample (see Table 1) thus consists of thirty-five ob- 
jects, almost equally divided between High Surface Brightness 
and Low Surface Brightness galaxies. 

For each galaxy we consider the exponential length h of 
the stellar disk and the radius Rq that we define as the radius 
at which the rotation curve reaches two thirds of its estimated 
asymptotic flat value Voo'. 



V(Rn) 



rVc 



(3) 



The velocity Voo is the estimate of the asymptotic velocity Voo 
obtained by averaging the rotation velocity in the outer parts of 
the rotation curve (see also Sect. 14.2b . The radius Rq thus mea- 
sures the slope of the inner approximately linear growth of the 
rotation curve. 



Name 


Type 


h/kpc 


i?n/kpc 


Refs. 


NGC 247 


HSB 


2.7 


3.5 + 0.6 


[4] 


NGC 300 


HSB 


1.8 


1.9 ±0.3 


[4] 


NGC 2403 


HSB 


2.1 


1.8 ±0.3 


[4] 


NGC 3109 


HSB 


1.35 


1.6 ±0.3 


[4] 


NGC 3769 


HSB 


1.5 


2.2 ±0.8 


[1], [2] 


NGC 3877 


HSB 


2.4 


2.6 ± 0.4 


[1], [2] 


NGC 3893 


HSB 


2.0 


1.7 ±0.4 


[1], [2] 


NGC 3972 


HSB 


1.6 


2.5 ±0.5 


[1], [2] 


NGC 3917 


LSB 


2.6 


2.6 ± 0.2 


[1], [2] 


NGC 4010 


LSB 


2.9 


3.8 ±0.8 


[1], [2] 


NGC 4085 


HSB 


1.3 


1.9 ±0.3 


[1], [2] 


NGC 4157 


HSB 


2.2 


2.4 ± 0.5 


[1], [2] 


NGC 4183 


LSB 


2.65 


2.9 ±0.8 


[1], [2] 


NGC 4217 


HSB 


2.4 


2.4 ± 0.5 


[1], [2] 


NGC 4236 


HSB 


2.7 


3 ±0.4 


[4] 


NGC 4389 


HSB 


1.2 


1.9 ±0.9 


[1], [2] 


UGC 128 


LSB 


6.8 


7 ±1.5 


[3], [5] 


UGC 1230 


LSB 


4.5 


6 ±1.2 


[3], [5] 


UGC 5005 


LSB 


4.4 


4.5 ± 0.6 


[3], [5] 


UGC 5750 


LSB 


5.6 


5.8 ±1 


[3], [5] 


UGC 5999 


LSB 


4.3 


5±1 


[3], [5] 


UGC 6399 


LSB 


2.0 


2.1 ±0.4 


[1], [2] 


UGC 6446 


LSB 


1.5 


1.6 ±0.6 


[1], [2] 


UGC 6667 


LSB 


2.45 


2.3 ±0.2 


[1], [2] 


UGC 6917 


LSB 


2.4 


2.3 ±0.2 


[1], [2] 


UGC 6923 


LSB 


1.1 


1.8 ±0.2 


[1], [2] 


UGC 6930 


LSB 


2.0 


2.4 ±0.3 


[1], [2] 


UGC 6962 


HSB 


1.4 


1.6 ±0.5 


[1], [2] 


UGC 6983 


LSB 


2.2 


1.9 ±0.5 


[1], [2] 


F 561-1 


LSB 


2.8 


2.8 ±0.5 


[3], [5] 


F 563-V2 


LSB 


2.1 


1.8 ±0.3 


[7], [6] 


F 565-V2 


LSB 


2.7 


3.5 ±0.5 


[3], [5] 


F 568-3 


LSB 


4.0 


3.1 ±0.5 


[3], [6] 


F 568-Vl 


LSB 


3.2 


2.7 ± 0.4 


[3], [6] 


F571-V1 


LSB 


3.2 


4±1 


[3], [5] 



Table 1. The final sample. The references are coded as: [1] Tully 
et al. (1996), [2] Verheijen & Sancisi (2001), [3] de Blok et al. 
(1995), [4] Kent (1987), [5] de Blok & McGaugh (1997), [6] 
Swaters et al. (2000), [7] McGaugh & Bothun (1994). 



E 
O 




Fig. 1. Correlation between the characteristic lengths h and Rn 
(see text for a description) for the sample of disk galaxies of 
Table 1 . The dashed lines represent 3cr error on the best-fit value 
of the angular coefficient. 
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When available (for example in the whole sample of galaxies 
belonging to the Ursa Major cluster), the reported exponential 
length is relative to the K-band photometry of the stellar disk, 
which is known to be the best representation of the actual stellar 
density distribution. In all the other cases we use the data as re- 
ported in the original source papers. We do not consider explicit 
error estimates for the exponential lengths h, since a uniform es- 
timate would not be available for all the galaxies in the sample. 
The estimate for the radius Rq is calculated directly, following 
its definition Eq. (O, from the observed rotation curve profile. 
The relevant error is estimated here by considering the errors in 
the velocity data-points and the process of interpolating the ro- 
tation curve between its measured points. 

Figure[T]demonstrates that there is indeed a strong, tight cor- 
relation between the two length scales considered. The correla- 
tion found is 



Rn/h^ 1.07 + 0.03. 



(4) 



The error on this measure is the standard statistical error, cal- 
culated by taking into account the uncertainties assigned to the 
values of Rq. 

Given the success of this test, it would be interesting to carry 
out a similar investigation on a better-defined, "cleaner" sample 
of galaxies, so as to test how much of the scatter around the 
correlation is intrinsic or rather due to the quality of the available 
data. However, a discussion of this point would bring us away 
from the main objective of this paper. 

The correlation coefficient for the entire sample is rh^R^ - 
0.96. It is thus interesting to note that the same correlation ap- 
plies to both HSB and LSB galaxies. The correlation provides 
additional support to a conspiracy argument, because it shows 
that even in the inner region of approximately linear growth of 
the rotation curve the dark matter distribution must be in tune 
with the distribution of the visible mass. Note that this is not 
necessarily in favor of the maximum-disk hypothesis, but rather 
it suggests that the relative distribution of visible and dark matter 
in spiral galaxies is homologous among different objects. Such 
homology could also be at the basis of observed correlations 
such as the Tully-Fisher (1977) relation (van Albada & Sancisi 
1986). 

In addition, this result will allow us to make a significant (al- 
though preliminary) test of the models developed in this paper 
on a simple idealized prototypical case that incorporates the ob- 
served correlation, without the need for separate tests on several 
observed objects. 



3. An isothermal halo embedding a zero-thickness 
disk 

We refer to standard cylindrical coordinates (R,z,0), so that 
r = yR^^+z^ represents the spherical radius. We assume that 
the rotation curve V{R) of the disk galaxy, for which we want to 
construct a mass model, is known (i.e., it has been "observed"); 
as a mathematical boundary condition, we assume that the rota- 
tion curve remains flat at Vca out to infinite radii. 

In order to make a fit to the observed rotation curve, we wish 
to construct axisymmetric equilibrium models made of two mass 
components: the visible disk and the dark halo. For simplicity, 
we imagine the disk to be made only of stars, with known surface 
brightness profile /i(^), and ignore the contribution of other visi- 
ble components. It will soon be clear that the contributions from 
a gaseous disk and from a spheroidal bulge with given structure 



could be easily incorporated, if desired. The visible disk is de- 
scribed as a zero-thickness disk; in a follow-up paper, we will 
address the study of a self-consistent finite-thickness disk as a 
generalization of the present model. 

The dark halo is taken to be characterized by a Maxwellian, 
isothermal distribution function. The choice of the dark matter 
distribution function is justified by an argument of physical sim- 
plicity. With the minimum number of free parameters, it guar- 
antees the flatness of the outer parts of the rotation curve. The 
procedure that we are going to describe can be easily adapted 
to study dark halos characterized by other physically justified 
distribution functions, if so desired (e.g., see Bertin, Saglia, and 
Stiavefli (1992) and Bertin (2000)). In conclusion, for the dark 
matter halo we take the isotropic distribution function: 



foMiE) 



Pdm 



(Ina-^) 



■ill 



expf-Zs/cr^j . 



(5) 



where £ = v^/2 + O^ is the specific energy of dark matter parti- 
cles and (i>j is the total gravitational potential. In our model Q>j 
is determined by the dark matter halo and by the disk: 



tr = ^DM + ^ 



^D ■ 



(6) 



The density distribution associated with the distribution function 
©is 

Pdm(R, z)^ fd\ foM = Pom exp [-Ot(R, z)/a^] , (7) 

where we have made it explicit that, in general, it is not spher- 
ically symmetric, because of the presence of the disk. At large 
radii, the gravitational contribution of the disk becomes vanish- 
ingly small, so that the total gravitational potential eventually 
approaches the natural spherical symmetry of the self-consistent 
isolated isothermal sphere. 

The dimensional free scales of the distribution function in 
Eq. (|5]) are the central density p^^ and the (constant) velocity 
dispersion cr. We identify the quantity p^^ with the central dark 
matter density, because we impose the condition 



Or(0,0) = 0. 



(8) 



In other words, we will only consider models that are regular at 
the center. 

The velocity dispersion cr is uniquely fixed by the value Vac 
of the rotation curve at large radii, so that the central dark mat- 
ter density p^^ is the only remaining free dimensional scale for 
the adopted dark matter density distribution. In fact, at large dis- 
tances from the center, where the system is almost spherically 
symmetric, matching with the observed rotation curve requires 



<ir(R,z)~<iT(r)~Viln 



ro 



which, inserted in the Poisson equation, gives 



y2 

' CO 

y2 



= 2(r^ 



AnGpl, 



(9) 



(10) 



We thus see that the scale length ro of the dark matter density 
distribution appearing in Eq. (|9]l is not a free length, because for 
assigned value of Voo (i.e., of cr) it is determined by the value of 
the central dark matter density, p^^- 
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For our zero-thickness disk, we take the surface density to be 
proportional to the observed surface brightness fJ.iR), i.e.: 



(M\ 
PD{R,z) = (-jlu(R)6iz); 



(11) 



here the function 6 is the Dirac delta-function. The ratio M/L is 
the mass-to-light ratio of the stellar disk, which we take to be 
constant; for assigned photometric profile fiiR), it represents the 
only free dimensional scale of the stellar density distribution. 

The model is constructed by solving the complete Poisson 
equation for the system, in which the matter densities are defined 
by the distributions Q and ( fTTT i: 



V^(bT{R,z) = 47tG [pd(R,z)+Pdm(R,z)] 



(12) 



This is a nonlinear differential equation for the total gravitational 
potential O^, which must be solved with the appropriate bound- 
ary conditions: Eq. ([8]l, as the condition at the center, and Eq. (|9]l, 
for the asymptotic behavior at large radii. The solution of the 
Poisson Eq. (fT2l i. which depends on the two free scales p^^ and 
M/L, will allow us to calculate the rotation curve of the model 
Vmod(R) - '<jR(d^T/dR)\,=o, to be compared with the observed 
V(R). 

Note that this method is intrinsically different from the stan- 
dard parametric decompositions of the mass models of spiral 
galaxies. In standard parametric analyses, the contribution of the 
dark matter halo to the total gravitational potential has a fixed 
shape (usually, spherically symmetric), which is independent of 
the presence and gravitational importance of the stellar disk and 
which can be changed only through a set of parameters a,, that 
is (^DM = ^DM{R,z;ai). The present self-consistent models do 
not have a fixed parametric profile for the shape of the gravita- 
tional potential of the dark matter component. Through Eq. (fT2l i. 
the dark matter gravitational potential, and thus its contribution 
to the total rotation curve, depends implicitly also on the prop- 
erties of the stellar disk (its surface brightness ^{R) and its mass 
to Hght ratio), that is O^m = ^dm{R, z;p%/^;fi(R), M/L). 

In conclusion, for a given rotation curve V{R) and surface 
brightness of the stellar disk //(/?), the best-fit self-consistent 
disk-halo decomposition will correspond to the self-consistent 
model defined by the pair (p^^^, M/L), together with the asymp- 
totic velocity Voo, that gives the best fit to the observed rotation 
curve. 



3.1. The mathematical problem in dimensionless form 

As units of length and of specific energy we will consider the 
scale length of the stellar disk h and the square of the velocity 
dispersion cr^ - y^/2, respectively, so that we introduce the 
dimensionless coordinates 



^ = R/h, ( = z/h,ri = r/h 

and the dimensionless gravitational potential 

•At- = ijJDM + ^D = -<^T/(yi/2) . 



(13) 



(14) 



We separate the Poisson Eq. ( fT2b in the two Poisson equations 
for the two different mass components of the system. The dimen- 
sionless equation for the stellar disk reads 



'}_d_ d_ 






lAo 



-^ t{^)5{0 



(15) 



Here we have introduced the normalized surface density, defined 
by t{R) = I.(R)/'L(0), and the parameter /3 



y8 = 



87rGS(0)/! M ^7rGfi(Q)h 



' or 






(16) 



which is the dimensionless version of the weight of the stellar 
disk as measured by the mass-to-light ratio M/L. 

To solve Eq. (fTSl ). we may use the standard technique of 
Hankel transforms (see Toomre (1963)). In this way the poten- 
tial associated with any surface density profile (M/L)fj.(R) can 
be calculated exactly. Consistent with the condition (jSj, we also 
impose that i/zoiO, 0) - 0. In particular, for purely exponential 
disks, £ = exp(-^), the associated gravitational potential is given 
by 



-A, 



..o-mj: 



d{hk) 



exp(-hk \0)Jo(hkO 



[l+(hkO^] 



3/2 



-1 



(17) 



where Jo follows the standard notation of Bessel functions. 

In turn, the Poisson equation for the dark matter halo can be 
written as 

Jgp^-gP + -gTij'f'DM = "« exp [(Adm(^, + <f'D(^,0] , (18) 



where the parameter a can be seen as the dimensionless expres- 
sion of the central dark matter density scale: 



' fVl 'A 



(19) 



The parameter a measures the halo concentration. With an eye 
to the central region of the system, a measure of the halo mass 
relative to the the mass of the disk is given by the parameter 
a/p — p^jjj^h/Y.0 (see also Eq. (fT2l) ). Thus the dimensional pair 
(p^j^,M/L) has its dimensionless counterpart in {a,/3). In the 
{a,jP) plane, the straight lines defined hy/3cca can be interpreted 
as lines of approximately given central density of the halo with 
respect to the mass of the disk; steep lines are those characteriz- 
ing relatively light halos. A better defined measure of the weight 
of the halo will be given in the next Section. 

The Poisson Eq. ( fTSl l for the dark matter halo has to be solved 
numerically (see Appendix A). Here the gravitational potential 
of the disk, calculated separately as described above, acts as an 
external potential that breaks the spherical symmetry and shapes 
the dark matter density distribution. The relevant boundary con- 
ditions can be obtained from the boundary conditions imposed 
on the total gravitational potential: 



^dm(0,Q) = 

^dm(ji) ~ -2 In (77 Vf ) - C ' '? » 1 



(20) 



where the symbol i//'^ indicates the constant value approached 
(with spherical symmetry) by the dimensionless gravitational 
potential of the disk at large distances from the center of the 

system. 
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4. Properties of self-consistent splieroidal 
isothermal halos 

The structure of the dark matter halo reflects the structure of the 
total gravitational potential, as prescribed by Eq. O. 

On the one hand, there is an artificial cusp in the vertical 
dark matter density distribution on the z -0 plane, which results 
from the singular (zero-thickness) density distribution assumed 
for the stellar disk. The detailed structure of the cusp can be 
obtained from the Jeans equation: 



Pdm{R,0) 



—pDM(R,0^)-—pDMiR,0-) 

oz oz 



d<l>D 



dz 



(«,()+) 



(S.0-) 



(21) 



This feature has no physical relevance and will be naturally elim- 
inated in the finite-thickness disk models that we plan to present 
in a separate paper. 

On the other hand, the halo is characterized by significant 
flattening. This latter feature is genuine and will not disappear in 
the regular finite-thickness case. Figure |2] shows the contours of 
the total gravitational potential, which coincide with those of the 
dark matter density, for the model defined by {a,l3) - (1 .9, 3.4), 
which, as we will see in the following, is the self-consistent 
model representative of the correlation in Eq. (|4]i. 

4.1. The shape of isothermal halos 

As naturally expected, the system tends to be spherically sym- 
metric in the outer parts (where the gravitational importance of 
the disk rapidly vanishes), but is highly oblate in the central re- 
gions. The flattening of the dark halo can be quantified in terms 
of the ellipticity e, defined as 






(22) 



where ^ and v(^) represent respectively the intersections with the 
^ = and ^ = planes of a given contour of the total gravita- 
tional potential. The profile of the ellipticity e as a function of 
radius is illustrated in Fig.|3] The high flattening inside one expo- 
nential length may be exaggerated because of the zero-thickness 
of the disk. Figure |4] shows the ellipticity reached by the self- 
consistent models at 2.2 exponential lengths in the {a,/3) plane. 
It would be interesting to compare the flattening obtained 
from our self-consistent models with the empirical determina- 
tions based on a variety of dynamical arguments, ranging from 
the dynamics of warped disks and polar rings to the flaring of gas 
layers (see the discussion in the paper by Oiling & Merrifield 
2000; see also the series of papers started by O'Brien et al. 
2010). It would also be interesting to compare the flattening of 
our self-consistent models with the indications that are expected 
to be gathered in the near future from weak and strong gravita- 
tional lenses (for example, see Trott et al. 2010); here, the natural 
comparison to be made is in terms of the relevant flattening pro- 
file of the projected density distribution. 

4.2. The flatness of self-consistent rotation curves 

A significant feature of the self-consistent models presented in 
this paper is the unexpected degree of flatness that character- 
izes their rotation curves. In principle, self-consistent rotation 
curves are designed to be flat only as an asymptotic condition 
for R » h, that is where the gravitational influence of the stel- 
lar disk becomes negligible. In the intermediate radial range the 
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Fig. 2. Contours of the dark halo density distribution (normal- 
ized to its central value) for the model defined by Eq. ( l30l l. Note 
that these contours coincide with those of the total gravitational 
potential. 



0.8 - 



0.6 - 



0.4 - 



0.2 - 




Fig. 3. Ellipticity profile of the dark matter halo illustrated in 
Fig. 12] as defined by Eq. (l22l l. 



shape of the rotation curve is instead determined by the mutual 
interactions between halo and visible disk and thus cannot be 
ascribed naively to the Maxwellian form of the halo distribution 
function (especially because in the intermediate radial range an 
isothermal sphere is known to be associated with a rotation curve 
exhibiting a characteristic oscillation pattern). 

To better demonstrate this point and to clarify the origin of 
the observed flatness, the self-consistent rotation curves will be 
compared with those of two different non-self-consistent mod- 
els: (i) the parametric models, defined by Eqs. ( 1112b ; (ii) the 
models, called for short semi-parametric, for which the field 
of the disk, with associated rotation curve defined by Eq. ([T]i, 
is superposed (non-self-consistently) to the field of a spheri- 
cal isothermal halo, i.e. a halo described by the density distri- 
bution in Eq. (|7]l, in which the total potential <i>TiR,z) is re- 
placed by the dark matter potential 'S>DMir). Therefore, any dif- 
ference in the rotation curves between semi-parametric models 
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Fig. 4. The ellipticity at 2.2 exponential lengths e(2.2) in the Fig. 5. Contours of the function (Eq. ( l23l l) for self-consistent 
{a,li) plane. models. 



and self-consistent models must be ascribed to the effects of self- 
consistency alone. 

With the aim of providing a systematic and quantitative anal- 
ysis, we start by defining as intermediate radial range the range 
2.2/1 < R < Ih, i.e. 2.2 < ^ < 7. Here the lower limit corresponds 
to the location at which the rotation curve of an isolated expo- 
nential disk reaches its maximum, while the upper limit repre- 
sents a reasonable choice of a location outside the bright optical 
disk. The quantity Vac will denote the average value of the rota- 
tion curve Vmod in the outer parts, 4 < ^ < 7. We then introduce 
the function 



©(«,/?) 



J2.2 



1 dV, 



mod 



Vo. d^ 



d^. 



(23) 



The function © is constructed so as to measure the flatness of 
the rotation curve V,„od in the intermediate radial range. The 
©-contours calculated from the self-consistent models are dis- 
played in Fig.|5] In Fig.|6]the two left panels show the contours 
of the function © obtained from the parametric models (upper 
frame) and the semi-parametric models (lower frame), for the 
same region of parameter space. A set of representative contours 
has been chosen arbitrarily, with the aim of describing the prop- 
erties of the three different cases. 

With respect to the range of values of the function ©, the 
parametric models and the semi-parametric models differ sig- 
nificantly from each other, because the semi-parametric models 
exhibit rotation curves that are less flat over the entire region 
of parameter space. This behavior is due to the fact that regular 
isothermal spheres are characterized by a well-known oscillation 
pattern in their density distribution and in the associated rotation 
curve, so that in the inner region the rotation curve overshoots 
its asymptotic value (see also Fig.|9]l. This property is removed 
a priori in the parametric profile of Eq. (|2|, often adopted in this 
context. In contrast, self-consistent models show a degree of flat- 
ness comparable to that of parametric models. The right panels 
of Fig. |6] illustrate (upper frame) the contours of the ratio be- 



tween the function © for parametric models and the same func- 
tion for self-consistent models and (lower frame) the contours 
of the ratio between the function © for semi -parametric models 
and the same function for self-consistent models. It is then clear 
that, in a relatively wide region of parameter space, where the 
contour levels exceed unity in the lower right frame of Fig. |6] 
self-consistent rotation curves are more flat than the coiTespond- 
ing semi-parametric rotation curves; as stated previously, this 
flattening effect can only be due to self-consistency (this point 
will be addressed further in Sect. 14.41 ). In addition, there is a re- 
gion of parameter space in which self-consistent rotation curves 
are more flat than parametric rotation curves (where the contour 
levels exceed unity in the upper right frame of Fig.|6l). This latter 
comparison has no special physical meaning, but it confirms the 
high degree of flatness of self-consistent rotation curves (which 
is also evident from Fig.|8]i. 

Another important point of comparison is the shape of the 
contours of the function © for the different cases. While the self- 
consistent models are flat in a well-defined region of parameter 
space (e.g., see the distinction between solid and dashed con- 
tours in Fig. |5]l, the parametric models are uniformly flat over 
an entire diagonal strip of parameter space (upper left frame of 
Fig.©. In Sect.|5]we will see that this strip is remarkably similar 
to the disk-halo degeneracy strip. 

4.3. The inner gradient of the self-consistent rotation curves 

After characterizing the properties of the self-consistent models 
at intermediate radii, in this subsection we focus on the proper- 
ties of the inner regions. We start by defining the radius Rn from 
the relation 



V(Rn) = 3^- . 



(24) 



with Vca introduced in Sect. 14.21 Figure |7] shows the contours of 
the quantity Rn/h for self-consistent models (black contours) 
and parametric models (gray contours). The thick contours in 



N.C. Amorisco and G. Berlin: Self-consistent nonspherical isothermal halos embedding zero-thickness disks 




Fig. 6. Top-left: Contours of the function © for parametric mass 
models. Top-right: Contours of the ratio between the function 
for parametric models and the same function for self-consistent 
models. Bottom-left: Contours of the function for semi- 
parametric models, with isothermal spherical halo. Bottom- 
right: Contours of the ratio between the function for semi- 
parametric models and the same function for self-consistent 
models. In the right frames, the thick contour represents unity. 
The axes in the four panels are {a,/3), as in Fig.|5] 



Fig.|2]represent the contours corresponding to the empirical cor- 
relation found in Eq. (|4|. 

The general shape of the contours is significantly differ- 
ent for the two models. While self-consistent models show 
roundish contours, parametric models exhibit a clear diagonal 
trend, which reinforces the expectation that a degeneracy prob- 
lem will arise in the corresponding disk-halo decompositions, as 
noted in the previous subsection. 

Finally, in the case of self-consistent models, the general 
shape of the Rn/h contours is similar to that of the contours 
of the function 0. The Rn/h value corresponding to maximum- 
flatness is slightly below the one suggested empirically by the 
analysis of Sect. 12.1! that is it corresponds to a slightly faster 
growing rotation curve (but see also the case of NGC 3198, de- 
scribed in Sect. 15.2b . 

4.4. Systematic properties of fhe self-consistent solutions in 
the available parameter space 

The aim of this subsection is to provide a description of the prop- 
erties of self-consistent rotation curves in the relevant parameter 
space. 

Figure [8] shows the rotation curves and the related disk- 
halo decompositions for an atlas of self-consistent models, in 
the region of parameter space of astrophysical interest (see also 
Sect.|5]l. In each frame the disk contribution to the rotation curve 
is dash-dotted, the contribution of the self-consistent halo is the 
solid black curve and the resulting rotation curve is indicated 
by a thick curve. Each frame records the pair {Mj,5,Mjx)). The 




Fig. 7. Contours of the quantity Rn/h, defined in Eq. (l24b . The 
gray contours, exhibiting a diagonal trend, are associated with 
the parametric models defined by Eqs. ([TJ and (|2|, while the 
black contours are associated with the self-consistent models. In 
both cases the thick line represents the empirical value found in 
Eq. (HJl, and the dashed lines represent the edges of the interval 
defined by a 3cr deviation. 



quantities M3 g and M7 indicate the ratio of the halo mass to the 
disk mass inside the spheres with radii of 3.5 and 7 exponential 
lengths, respectively. 

As already noted in Sect. 14.21 in general the rotation curves 
are significantly flat and lack specific features. The flattening ef- 
fect described in Sect. 14.21 shows up clearly in Fig. [8] in a sort 
of suppression effect. Here the effects on the self-consistent halo 
caused by the presence of the stellar disk can be followed by 
looking at each individual column, that is by fixing the value of 
the dimensionless central density of the halo a and by letting the 
weight of the disk /3 increase. The parameter a fixes the inner- 
most gradient of the dark halo contribution to the rotation curve 

(which is proportional to A/Gp^^). However, in general self- 
consistency cannot support a fast growing contribution from the 
dark halo in the presence of an important stellar disk, and then 
it "softens" it near 1 H- 2 exponential lengths. At small values 
of /3 self-consistency suppresses the bump in the halo contribu- 
tion to the rotation curve characteristic of the isothermal sphere 
(e.g., compare the (a = 9,y6 = 0.5) frame to the {a = 9,/3 - 2) 
frame). For heavier disks (large values of /3), self-consistency 
suppresses the halo contribution at intermediate radii (e.g., com- 
pare the (a = 0.5, /3 - 3.5) frame to frames corresponding to 
higher values of /3). For high values of both a and /? (see for 
example {a - 9,/3 = 8)) self-consistency makes the halo con- 
tribution fairly constant, but with a value which is significantly 
below 100 %, which will be reached eventually, at much larger 
radii. 

This suppression effect has three main consequences on 
the general shape of self-consistent rotation curves: it flattens 
them globally (with respect to the rotation curves of the semi- 
parametric models; see Sect. I4.2l i: it removes the oscillations 
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characteristic of the isolated isothermal sphere (provided the 
stellar disk is not too light); it makes the signature of the pres- 
ence of an important luminous disk more evident, by lowering 
the halo contribution. 

Similar general conclusions can be made by inspection of 
Fig. |9] Total rotation curves of self-consistent models (thick 
gray lines), parametric models (black solid lines), and semi- 
parametric models (dashed lines) are shown here in comparison, 
in a general format and layout similar to those of Fig. [8] In each 
frame, the stellar disk contribution is given by the dash-dotted 
line, while the constant line V„,od - ^oo is drawn to make the 
comparison easier. What clearly catches the eye is the remark- 
able uniformity of the parametric rotation curves, when com- 
pared to those of the semi-parametric models and of the self- 
consistent models. From a physical point of view, this unifor- 
mity is not justified, because it confirms that the parameteriza- 
tion used flattens out arbitrarily the properties of the physical 
isothermal solutions that it is meant to represent approximately. 



occurs in the vicinity of a single point of parameter space. We 
are then brought to conclude that degeneracy may be broken by 
self-consistency. This last point is confirmed and better quanti- 
fied in the following separate Section. 



4.5. Conspiracy and degeneracy 

Figure|7]demonstrates that the available parameter space cannot 
accommodate values of Rn/h significantly different from unity, 
especially if we also require that the rotation curve be flat at in- 
termediate radii. If galaxies were found with Rq too small or too 
large with respect to h, there would be no way for the picture 
of disks embedded in quasi-isothermal halos to hold and then 
one should be led to look for dark matter distributions signif- 
icantly different from those characteristic of isothermal halos. 
This point will be further addressed in Appendix B, where the 
difficulty of fitting such "anomalous" rotation curve is demon- 
strated. Fortunately, the observed correlation Eq. (|4|l shows that 
the picture of a quasi-isothermal halo is indeed viable. 

Figure [8] confirms that self-consistent models lead to rota- 
tion curves that are realistic, in the sense that they are smooth 
and featureless for a wide range of disk-to-halo mass ratios; 
note that the maximum disk solution associated with the cor- 
relation (HJi (see Sect. l5.1l for description) is expected to occur at 
/3 ~ 5 (cf. Eq. (l28Tl), approximately in the middle of the param- 
eter range discussed and explored in this Section. An inspection 
of Fig. |9] shows that parametric models are also characterized 
by this property, actually in an even more prominent way. In a 
sense, these findings demonstrate that we are not facing an issue 
of conspiracy or, at least, that the conspiracy is not as dramatic as 
often stated. In fact, differently from what is naturally expected 
and is often reported, no fine tuning is actually required in this 
problem, because for very different values of the halo-to-disk 
mass ratio the rotation curves resulting from the combination of 
a disk with a quasi-isothermal halo, turn out to be (surprisingly) 
smooth and flat. In this respect, the correlation represented by 
Eq. (4) implies a much tighter coupling between dark and lumi- 
nous components in the {a,/3) plane. 

In the context of parametric models, these conclusions ap- 
pear to be quite artificial. In fact, the use of the simple model 
associated with Eq. ^ removes a priori significant features that 
are characteristic of the regular isothermal sphere. Such non- 
natural uniformity in parameter space ends up in generating the 
problem of disk-halo degeneracy. In fact, the diagonal trend of 
the parametric models exhibited in the upper-left frame of Fig.|6] 
and in Fig.|2]reinforces the diagonal trend exhibited by the same 
models in Fig. [TOl thus anticipating a problem of degeneracy of 
solutions along diagonal lines in parameter space. 

The self-consistent models, instead, have a distinctly differ- 
ent behavior (see Fig. |5] and Fig.]?); for them, maximal flatness 
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5. Self-consistent disk-halo decomposition: the 
removal of the disk-halo degeneracy 

5.1. A preliminary test on a fiducial rotation curve 

As anticipated briefly in Sect. 12.11 the empirical correlation (|4]i 
found between the two relevant scale lengths h and Rn allows us 
to make a preliminary test by focusing on a single, representa- 
tive idealized case. Therefore, a discussion of specific observed 
objects, such as NGC 3198, is postponed and will be carried out 
separately (see following Sect. I5.2l i 

As a typical reference case, we consider a galaxy character- 
ized by a purely exponential stellar disk, 1,(R) = exp(-^), and by 
the following rotation curve 



y(^,T)/ye„ = i-exp(-^/T) 

If we take 

T~ 1 , 



(25) 



(26) 



the correlation found in Eq. (|4|l is reproduced. 

The simple analytical form of the adopted rotation curve 
allows us to calculate explicitly the value of the "observed" pa 
rameter Oobs defined in Eq. (IA.13b : 



aobs 



aobsir) = exp(2y)/2T 



(27) 



where y s; 0.577 denotes the Euler gamma constant, thus 
o^ohsiT = 1) ~ 1.59; the following discussion will be based on 
the use of this value of agbs- 

Note also that, for a given value of t, the rotation curve 
( l25T l will admit a well-defined maximum-disk decomposition. 
In practice, before carrying out any detailed decomposition, we 
may introduce, as a simple definition of maximum-disk, the disk 
characterized by the value of the dimensionless weight Pmaxi^) 
that gives, without any dark matter contribution, the best fit to 
the inner rotation curve in the interval [0,/?n]- Note that this 
maximum-disk decomposition is designed in such a way to re- 
cover the central gradient of the rotation curve. Sometimes a dif- 
ferent definition is used, by referring to the stellar disk which 
reaches the asymptotic rotation velocity in correspondence of 
the maximum of its rotation curve. This definition gives higher 
values of yS {Pnmx ~ 10 for an exponential disk), but it is com- 
pletely independent of the properties of the central part of the 
rotation curve. 

For T given by equation ( l26l l. we find 



A»«.(l) ^ 5 ; 



(28) 



this value of j0 identifies an important reference value to which 
the results obtained from the self-consistent decomposition will 
thus be compared. We recall that /? = 5 is approximately in 
the middle of the parameter range discussed and explored in 
Sect. 1441 

The properties of the self-consistent models constructed in 
the present paper are best illustrated by describing their ability 
to fit the just described representative idealized case in compar- 
ison to a fit performed by more standard parametric analyses. 
For both methods the goodness of a disk-halo decomposition, 
defined by a pair ia,/3), is quantified by the function 



H(a,yS,T) 



^1^"^-"' 




Fig. 10. Contours of the function H for the case of the parametric 
decomposition based on Eqs. ^^. As in similar plots (Figs.fTTl 
IB.1IIB.2IIB.4IIB.5I ). the displayed contours are chosen arbitrarily 
to illustrate the properties of the H function. 



that is the integrated squared residuals between the "observed" 
rotation curve and the rotation curve calculated from the model 



(^)-V(^,T)f d^.. 



(29) 



V„,„d - ■JV^ + V^^. The cut of the integration at seven expo- 
nential lengths is a reasonable choice, consistent with the anal- 
ysis of Sect. 14.21 . In this preliminary test we do not treat the 
asymptotic velocity as an additional free parameter of the model. 
In other words, we just analyze the deviations defined by Eq. 
(29) in the (a,/3) plane, at fixed Vca- The asymptotic velocity 
will be kept as a free parameter in Sect. 5.2, for the case of NGC 
3198. 

The results obtained from the standard parametric decompo- 
sition are shown in Fig. [10] which clearly exhibits the disk-halo 
degeneracy pattern. In particular, consider the contour marked 
by 0.02. A sizable diagonal strip in parameter space consists of 
points that are basically equivalent from the point of view of the 
quality of the fit, even though they correspond to models that 
are physically very different, ranging from fairly light disks up 
to maximum-disk solutions (with /3 changing by a factor higher 
than ten, from « 0.5 to x 6.5). 

In contrast. Fig. [TT] shows the contours of function ( |29] l for 
the case in which the model rotation curve is calculated using the 
self-consistent method presented in this paper The plot ranges 
and the values of the solid line contours plotted in Fig.fTOland 
Fig. [TT] are the same. The self-consistent models give a better 
disk-halo decomposition for three distinct reasons: firstly, the 
value Ebf = E(abf,Pbf, 1 -02) ^ 0.00056 of the function H corre- 
sponding to the best-fit model (abf,/5bf) is more than ten times 
smaller than the value attained by H based on the simple para- 
metric decomposition (HJ'!!^ » 0.006); secondly, the area of the 
parameter space inside the contour 2S/,/ » 0.001 12 (or any other 
sensible multiple of H/,y) is quantitatively smaller, that is the 
best-fit model is identified more sharply; finally, the contours 
of the function H in the self-consistent case do not show any di- 
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Fig. 11. Contours of the function S for the case of the self- 
consistent decomposition presented in this paper. 
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Fig. 12. Disk-halo decomposition associated with the best self- 
consistent fit. The thick dashed line is the "observed" rotation 
curve, the thick solid line is the best-fit rotation curve, the thin 
lines represent the disk (dash-dotted) and the halo (solid) contri- 
butions. 



agonal trend, thus confirming that the disk-halo degeneracy is 
removed. 

In conclusion, the best-fit model identified by the self- 
consistent analysis is defined by the following values of the rel- 
evant dimensionless parameters 



(afc/,/3i/)^ (1.9,3.4). 



(30) 



Therefore, the best-fit model is significantly sub-maximal (see 
Eq. (l28ll). The associated rotation curve is shown in Fig. [12] 

In order to better demonstrate the beneficial role of self- 
consistency in breaking the disk-halo degeneracy, we now con- 
sider the self-consistent models that coiTespond to parameters 
with good (degenerate) parametric fits to the "observed" rota- 
tion curve, but are ruled out by the self-consistent analysis. In 




Fig. 13. Disk-halo decomposition by the (ahm^Phm) self- 
consistent model. The format and the notation are the same as 
in Fig. [12] The rotation curve of the corresponding {akmiPhm) 
parametric model is displayed as a thick dashed gray line. 



V 




Fig. 14. Disk-halo decomposition by the {aim,Pim) self- 
consistent model. The format and the notation are the same as in 
Fig. [12] The rotation curve of the coiTesponding (a/m, /?;,„) para- 
metric model is displayed as a thick dashed gray line. 



particular, as representative cases we choose the high-mass (i.e., 
high-y6) and low-mass (i.e., low-yS) ends of the 2HJJl'^ contour in 



(ahnuPhm) 
(aim.Plm) 



(0.9,5.4) 
(6,1.75) 



(31) 



The self-consistent rotation curves associated with these two 
models are shown in Figs. [13] and [14] Note that the diff'er- 
ences between the calculated total rotation curve and the "ob- 
served" one are significantly larger than any realistic eiTor-bar. 
Furthermore, the rotation curve of the high-mass case confirms 
the properties of the self-consistent rotation curves described in 
Sect. 133] 



5.2. Application to ttie case of NGC 3198 

As a specific example, we now consider the classical case of 
NGC 3198. We refer to the rotation curve reported in Table 2 of 
van Albada et al. (1985) and to their measure of the exponential 
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length (h - 2.68 kpc). The value of the ratio Rn/h estimated 
from the rotation curve is Rn/h - 0.6 ± 0.1. This is on the low 
side of the correlation reported in Sect. 12.11 that is the rotation 
curve of NGC 3198 rises more rapidly than average in its cen- 
tral parts. This feature may be partly related to the fact that the 
photometric profile of NGC 3198 is not exactly exponential in 
the innermost parts of the disk (see Begeman 1989; see also the 
discussion by Dutton et al. 2005). For simplicity, in the follow- 
ing we will consider the galaxy as bulgeless and we will ignore 
the contribution of the gaseous disk. We recall that the simple 
analysis provided in Sect. 15.11 applied to the inner parts of the 
rotation curve, suggests that the maximum-disk decomposition 
for NGC 3198 has/?;^,^c3i98 ^ ^q^ 

In this case the rotation curve is provided with its observa- 
tional errors. To evaluate quantitatively the goodness of a given 
disk-halo decomposition, we then adopt the classical likelihood 
function: 



A(Va.,a,j3)^expL^Yj 



^^[V,„od(Rd-V(Ri)f 



Var[y(/?,)] 



(32) 



where Var [V(7?,)] denotes the square of the error associated with 
the /-th rotation curve data-point. In order to bring out the ex- 
plicit dependence of the likelihood A on the assumed value of 
Voo, we note that the previous equation can be written as: 



^ /=28 



[Va.-v„oARi,a,/3)-V{Rdf 
Wai[V(Rd] 



,(33) 



in which the function v,nod{R, <^,P) is normalized to unity at large 
radii. 

We note that the rotation curve data from four to eleven ex- 
ponential lengths suggest 



Voo^ (149.8 + 0.5) km/s. 



(34) 



However, in order to avoid a possible bias that might be intro- 
duced by "freezing" this quantity, we treat the asymptotic ve- 
locity Vco as a free parameter, to be determined by the fitting 
procedure together with a and /?. 

Figure [15] shows the results of the parametric disk-halo de- 
composition in the {a,l3) plane, that is, after integrating out the 
dependence on Voo in the likelihood Eq. (l32l l: 



A(a,y6) 



X" 



fi?VocA(Vco,Q',yS). 



(35) 



The contours identify regions of 68% and 95% confidence level; 
the best-fit model is marked by a full dot. Interestingly, two max- 
ima exist, located along the disk-halo degeneracy strip described 
in the previous subsection, but with significantly diff'erent phys- 
ical properties. In fact, for varying values of the asymptotic ve- 
locity Voo, the location of the best-fit model in the (ff,/?) plane (at 
fixed Voo) changes. As a result, the integrated likelihood Eq. (l35T l 
actually increases the amount of degeneracy associated with the 
parametric disk-halo decomposition. This is particularly evident 
from Fig. [16] which shows the projections along the two sepa- 
rate parameters a and /? of the the likelihood itself, exhibiting a 
clear bimodality. The disk-halo decomposition is unable to se- 
lect between a solution very close to the maximum-disk and a 
decomposition with a much lighter disk. The coiTesponding de- 
compositions are illustrated in Figs. [17] and [18] which show the 
two best-fit parametric decompositions, with projections on the 
(a,y8) plane contained respectively in the two separate areas of 



y0 




10 



15 



a 



Fig. 15. Contours of 68% and 95% confidence regions in the 
(a,l3) plane according to the likelihood defined in Eq. ( 1351 ) for 
the parametric decomposition. The full dot marks the best-fit 
model. 



the confidence regions in Fig. [15] This curious result may be re- 
lated to the specific structure of the rotation curve of NGC 3 198 
(but see also Appendix B). In any case, the general statement 
holds, that the parametric decomposition is unable to lead to a 
unique determination of the weight of the stellar disk. 

For completeness, we record the coordinates of the best-fit 
model in the (Voo, a, (3) space, together with the related reduced 



'^ value: 



(145.9,8.7,3.85,0.6). 



Here we have defined 



2 _ 
X red 



1 



28-3 



!=28 

Z 



Wai-[ViRd] 



(36) 



(37) 



The situation is completely different if the disk-halo decom- 
position is performed by means of the self-consistent models 
presented in this paper Figure[T9]shows the contours of the 68% 
and 95% confidence levels and the best-fit model, calculated us- 
ing the same likelihood as for the parametric case; note the sig- 
nificant change in their size compared to Fig. [15] At fixed values 
of the asymptotic velocity Vca , the general behavior of the best-fit 
models is similar to the one noted for the parametric analysis: the 
higher the adopted value of Vca, the higher the value of/3 found in 
the best-fit decomposition. However, here the differences cancel 
out in the integration Eq. (l35b . because of the different sharpness 
of the best-fits obtained at fixed values of Vca- As a result, the 
self-consistent decomposition is determined uniquely. In Fig. [20] 
the likelihood projections are indeed smooth and do not exhibit a 
bimodal shape. As for the parametric decomposition, we list the 
coordinates of the best-fit model in the (Vco, a,P) space, together 
with the related reduced ;\f^ value: 



(160.1,2.6,5.9,0.7). 



(38) 
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Fig. 16. Projections of the likelihood (l32l i for the single parame- 
ters (a,y6) for the parametric decomposition. 
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Fig. 17. Parametric disk-halo decomposition of NGC 3198 cor- 
responding to the maximum of the likelihood in Fig. [15] with 
higher disk importance. The coding is the same as in Fig. [12] 



V 
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Fig, 18. Parametric disk-halo decomposition of NGC 3198 cor- 
responding to the maximum of the likelihood in Fig. [15] with 
lower disk importance. The coding is the same as in Fig. [12] 
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15 



Fig. 19. Contours of 68% and 95% confidence regions in the 
{a,l3) plane according to the likelihood defined in Eq. ( l32b for 
the self-consistent decomposition. The full dot marks the best-fit 
model. 




Fig. 20. Projections of the likelihood ( l32l l for the single parame- 
ters (a,/?) for the self-consistent decomposition. 



The best-fit value of the asymptotic velocity Voo is higher than 
the value suggested by the rotation curve in Eq. (l34l l. consistent 
with the properties of the rotation curves displayed in Fig.|8]and 
with the suppression effect described in Sect. 14.41 

We note that the reduced x' value of the best self-consistent 
fits is somewhat higher than for the corresponding parametric 
fits. However, this difference is small and the reported values 
for the reduced x^ account for globally good fits. This is evi- 
dent from Fig.|2T] which shows the best-fit self-consistent disk- 
halo decomposition together with the observed rotation curve. 
However, the model has a certain difficulty in fitting the second 
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Fig. 21. Self-consistent disk-halo decomposition of NGC 3198. 
The coding is the same as in Fig. [121 



and third data-points, which remain higher than the best-fit self- 
consistent rotation curve. Clearly, this could be ascribed to the 
exponential modeling of the density profile of the visible matter, 
which, as noticed before, does not take into account the small 
bulge and the gas present in NGC 3198. Since these two data- 
points are those that define the value of the ratio Rnjh for NGC 
3198, Fig. Q confirms that the best-fit self-consistent model has 
a slightly higher value for this ratio. Note also that the proper- 
ties of the best-fit self-consistent decomposition are quite differ- 
ent from those found in the parametric analysis, with the best-fit 
model characterized by a disk of intermediate weight. As is ev- 
ident from Fig. |2T] the disk contribution to the rotation curve 
is important in the inner parts of the system. However, with its 
jS X! 6, the best-fit self-consistent model remains below the clas- 
sical maximum-disk solutions that van Albada et al. (1985), for 
example, place at/? ^ 9. 



6. Conclusions 

This paper has addressed the construction of self-consistent 
models for a dark matter distribution assumed to be isother- 
mal, in the presence of a zero-thickness disk. The models have 
then been applied to interpret the properties of observed rotation 
curves, with special attention to the issues of conspiracy and de- 
generacy that are known to affect their disk-halo decomposition. 

We argue that what has been learned from the study of the 
strictly isothermal case should be applicable to the case of dark 
matter distributions characterized by a quasi-isothermal inner re- 
gion (as is the case for the results of coUisionless violent relax- 
ation). 

It would be interesting to test the results on the dark matter 
distributions predicted by current cosmological models, which 
definitely favor a density profile asymptotic to r^^ . This test 
could be performed with an approach similar to the one adopted 
in this paper, provided that we had a physically justified distribu- 
tion function to start from. In particular, it would be interesting 
to test to what extent the cosmologically predicted distributions 
could be reconciled with the existence of cases, such as that of 
NGC 3198, in which the rotation curve remains very flat in a 
wide radial range. 

Note that this paper has considered the impact of a given visi- 
ble matter distribution on the distribution of dark matter A study 
of the dynamical role of the dark halo in shaping the distribution 



of the visible matter would require considering a formation sce- 
nario that we do not wish to address here. 

The main results obtained in this paper are the following: 

- As a preliminary empirical study, we have shown that for 
both high surface brightness and low surface brightness disk 
galaxies the scale length h of the light profile of pure disks 
correlates with the scale length Rq (defined in Eq. (O) that 
characterizes the central growth of the rotation curve. This 
suggests that the distribution of dark and visible matter in 
pure disk galaxies is homologous, i.e. characterized by simi- 
lar global relations in different galaxies, even if significantly 
different in size and type (HSB vs. LSB galaxies). Based on 
the results obtained in the rest of the paper, this point is inter- 
preted as an indication of a gravitational interplay between 
luminous and dark components in the innermost parts of the 
system. We argue that this homology is more relevant and 
surprising as a conspiracy than the actual flatness of the ob- 
served rotation curves, because, especially in the context of 
parametric models, rotation curves turn out to be almost flat 
for very different values of the disk-to-halo mass ratio. 

- We have constructed self-consistent models of isothermal 
halos embedding zero-thickness stellar disks, thus obtain- 
ing fully self-consistent axisymmetric isothermal dark mat- 
ter halos. These make basically a two-parameter family of 
models; the two parameters measure the weight of the disk 
and the central density of the halo (see Eqs. ( fTglfTSl l). A third 
parameter, the asymptotic value of the rotation curve at large 
radii, is in principle free but can often be treated as a fixed 
parameter, given the fact that it is generally well constrained 
by the observations; note that a discussion with the third pa- 
rameter kept free is given in our study of NGC 3198 (see 
fifth item in this subsection). The adopted solution procedure 
can easily incorporate the presence of a zero-thickness disk 
of gas (with desired surface density). In addition, a similar 
procedure could be followed to construct models based on 
different choices of the dark matter distribution function. In 
turn, we can also argue that the qualitative conclusions ob- 
tained in this paper by calculating the properties of strictly 
isothermal halos should also apply to more general systems 
where the halo distribution, in the region of interaction with 
the disk, is quasi-isothermal. 

- We have shown and quantified that, for astrophysically re- 
alistic values of the two free parameters of the models, the 
dark matter halo is significantly flattened, especially in the 
regions where it coexists with the bright optical disk. Inside 
one exponential scale length, the dark matter distribution is 
that of a thick disk; nevertheless, its contribution to the rota- 
tion curve is approximately linear in the same region, denot- 
ing the presence of a natural core. 

- The rotation curves of the self-consistent models have been 
shown to be generally smooth and flat in their intermedi- 
ate regions. In particular, self-consistency has been proved 
to provide a physical "mechanism" able to flatten rotation 
curves in the region of the visible disk. This flattening effect 
is stronger in the region of parameter space that appears to 
represent realistic systems. This result goes against the naive 
expectation that self-consistency should tend to enhance fea- 
tures present in the potential well and then in the related ro- 
tation curve. 

- We have shown that the application of self-consistent mod- 
els to decompose an observed rotation curve breaks the disk- 
halo degeneracy. This result has been achieved under some 
simplifying assumptions. First, for simplicity we assumed 
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that the density profile of the disk is a perfect exponential 
(although our self-consistent models can be constructed for 
any zero-thickness axysimmetric density profile for the stel- 
lar disk). Second the dark matter distribution function has 
been assumed to be a perfect isothermal; this point can be re- 
leased in favor of other physically justified choices. The disk- 
halo degeneracy breaking has been obtained first by studying 
an idealized rotation curve representative of a typical rota- 
tion curve of a pure disk spiral galaxy and then by a direct 
application to the classical case of NGC 3198. We have also 
provided arguments that support the idea that, to some ex- 
tent, such degeneracy problem has been artificially created 
by the use of parametric models (such as the model associ- 
ated with Eq. (|2])). In fact, the rotation curves generated by 
these parametric models have been shown to be exceedingly 
uniform and featureless. 

- We have thus derived the best-fit self-consistent disk-halo 
decompositions and proved that they favor mass decompo- 
sitions characterized by stellar disks of intermediate weight, 
which are still important in the central regions of the sys- 
tem, but are significantly lighter than maximum-disk solu- 
tions. On the other hand, the homology property noted in the 
first item of these Conclusions would still make such sub- 
maximal solutions consistent with the small scatter of the 
Tully-Fisher (1977) relation. Of course, a full demonstration 
of this conclusion would require that we moved beyond the 
case of the fiducial rotation curve (see Sect. 5.1), by study- 
ing, in the way described in Sect. 5.2 for NGC 3198, a com- 
plete sample of spiral galaxies. But this would obviously go 
well beyond the objectives of the present paper 

- The resulting rotation curve has been interpreted as resulting 
from a suppression effect, characteristic of self-consistency. 
Self-consistency would affect a proposed maximum-disk de- 
composition by significantly reducing the dark halo contri- 
bution to the total rotation curve and thus leading to worse 
fits, especially at intermediate radii. 

- The resulting M/L values for the stellar disk should be stud- 
ied for each individual galaxy, in view of their implications 
and consistency with respect to what is known about the rel- 
evant stellar populations. In this context, an important test 
would be the study of LSB galaxies. However, this very im- 
portant issue goes well beyond the goals of the present paper 

In a follow-up paper we will generalize the results of this 
paper to the case in which the disk is characterized by finite 
thickness, with results that may be of interest to projects such 
as the Disk Mass Project (Verheijen et al. 2004 and Westfall et 
al. 2008). 

Acknowledgements. We wish to thank Luca Ciotti, Wyn Evans, Renzo Sancisi, 
Tjeerd van Albada and an anonymous Referee for making a number of interest- 
ing comments and suggestions. 

Appendix A: Construction of the models by an 
iterative procedure 

Equation ( fTSb . with the boundary conditions ( |20] |. is solved by 
means of a modified version of the iterative procedure presented 
earlier by Prendergast & Tomer (1970) in view of a completely 
different problem (the study of nonspherical models for rotat- 
ing elliptical galaxies). It is based on the classical expansion 
in Legendre polynomials, a natural approach in the case of ax- 
isymmetric systems. However, the original method is viable only 
for gravitational potentials that are regular at large radii. In con- 
trast, in our model, the potential of the isothermal halo diverges 



at large radii. This difficulty has been overcome by separating 
the potential of the dark matter halo into two different parts: a 
regular part, which converges to finite values at large radii, and 
a divergent one, which follows the asymptotic prescription of 
Eq. ( l20t . Below we provide the details of the technique devel- 
oped in the present paper. Here it is convenient to work with 
standard spherical coordinates (r, <p, 6) and their dimensionless 
analogue (rj, </>, 6). We start by recalling briefly the method de- 
vised by Prendergast & Tomer (1970). The solution of the two- 
dimensional Poisson equation 

VV('7,0)=P['7,0;<A(^,0)] (A.l) 

is obtained by iteratively solving the equations 

VV*"^"('7, 0) = P [ri, 0; ip^"\ri, 61)] (A.2) 

to which the same boundary conditions as in Eq. dA.l! ) are im- 
posed. The potential at the iterative step (n -i- 1), i/r*"+'\ is ob- 
tained from the potential at the previous step, \l/^"\ by solving 
Eq. ( IA.2! ) exactly, using the standard multipole expansion in 
Legendre polynomials, here denoted by Pk'- 



p("'(/7,0) = 2pf(7;)A(cos0) 



k=Q 



and 



^<"+'>(77,0) = £ i^^f(rf}Pu(co^e) 



(A.3) 



(A.4) 



A=0 



For completeness, we display the general solution of the «-th 
iterative step: 

r]'j^^r]"-'pf\rj')di]'+ (A.5) 

f(l]')dTj'] 



Y'^'o P^ (cos ^) 

~ 2ji:=2 2k+l 

^if*' Jo 'I f^k 



where we have used the notation i^(0, 0) - T. The iteration can 
be stopped when the desired accuracy prescription is met, as for 
example when in the relevant domain D we find that 



max 



(n+l) _ ^(") 



^(n+1) + ^(«) 



< £ . 



(A.6) 



The dark matter gravitational potential is then separated into two 
parts: 

I^DM = Ipasy + 4f ■ (A.7) 

The potential defined as if^asy obeys the Poisson equation 

VVflJJ = Pasy (A. 8) 

and has the same asymptotic behavior as in Eq. ( l20b . In this way 
the potential tfr converges to zero at large radii and the related 
Poisson equation 

7 g^^gZ + ^ U = "'^ ^^P [^"^y + ^D + 4']- Pasy (A.9) 

can be directly solved by the iterative multipole expansion out- 
lined previously. 

Obviously, there is an infinite number of pairs {pasy, 4'asy) 
that meet the conditions of Eqs. ( IA.8b and ( l20b . Therefore, the 
construction method proposed below just reflects one reasonable 
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choice. We decided to take pasy and tpasy with spherical symme- 
try, because the relevant asymptotic condition for ij/asv is charac- 
terized by spherical symmetry and because the spherical Poisson 
equation admits a simple explicit solution: 



Jo ^ Jo 



f^asyirf) = lAfl.n'(0) + 



The precise choice of the form of the density pasy is guided by the 
goal of simplifying the following numerical procedure to solve 
Eq. ( IA.9b . Different choices of the density poj,, correspond to dif- 
ferent shapes of the potential for the same pair of free param- 
eters (a,/3). The adopted practical recipe to construct the density 
Pasy turns out to be useful and efficient. 

We start by defining what we call an "observed" pseudo- 
potential, constructed from the observed rotation curve: 



<i>obs(r) 



f 

Jo 



VHs) 



ds . 



At large radii it has the asymptotic expression 



(A. 11) 



(A.12) 



where the radius Tq*' can be calculated from the precise form of 
the rotation curve. The behavior of Eq. ( IA.12I ) naturally relates 
the potential (b„bs (and its dimensionless counterpart ijjohs) to a 
density distribution with isothermal functional shape: 



\ d 2 d 
Tj'^ dr] drj 



-a„bs exp [il/obsir])] , rj » I 



(A.13) 



As shown in Eq. ( fTOt . the value of the constant aghs is deter- 
mined by the value of the radius rg*' only, and is thus fixed. 
Having made this points clear, we can define /)„„. to be 



Pasy(j]) = -a„hs exp 



^obs 



s[^ff^ 



+ P^°'W: 



(A. 14) 



where we have introduced a regularizing core structure with 
characteristic size tJc and subtracted the monopole term of the 
stellar density 



*"4X! 



I;(77) 6{o A)(cos e)dcose = - — 

2 T] 



(A. 15) 



With this density profile, from Eq. ( I A. 10b we calculate 
the gravitational potential tpasy, using the free constant i/'n.sy(0) 
to make it meet exactly the asymptotic behavior required by 
Eq. ( l20b . The determination of the correct value of the constant 
t//asy(0) determines also the boundary condition to be used for 
Eq. ( IA.9b . which is set by the prescription: 

MO, 0) = iA(0, 0) + i/,asy(0) + iAd(0, 0) = . (A. 16) 



A.1. The truncation of the angular expansion 



A code has been written to solve Eq. (IA.9I) through the technique 
described in the previous section. The code is able to manage au- 
tomatically the appropriate number of multipole terms to meet a 
required accuracy (see Eq. ( IA.6b ). At each iterative step, the rele- 
vant number of multipole orders is calculated from the following 
prescription: 



iAf(>7)P^(cosg) 



— > 10 



i-i. 



(A. 17) 




too 



Fig.B.l. Contours of the function H for the parametric decom- 
position (based on Eqs. ([T] |2]l) of the rotation curve defined by 
Eq. dB]) with Tf = 0.4. 



In words, at the «-th iteration step, the multipole term of order 
k is retained in the expansion if its (relative) contribution to the 
gravitational potential (/r*"' (calculated with the previous k - I 
multipole orders only) is comparable to the accuracy level e we 
want to meet. This comparison is made at the coordinates (^, 6), 
chosen in the region where the deviation from spherical symme- 
try is strongest: x n/2 and 77 ^ 1 ^ 2. 

The models described in the paper are calculated with an 
accuracy prescription of e = 10""*. Obviously, the number of it- 
erations required by such precision depends on the gravitational 
importance of the stellar disk (which determines the flattening 
of the halo). For values of the pair ia,p) that represent an as- 
trophysically realistic configuration (see Sect. |4]l, the number of 
required iterations is 10 H- 20, with a number of multipole orders 
of 12 H- 20. 



Appendix B: Anomalous rotation curves 

In this Section we will analyze the case of rotation curves 
with inner gradients that are significantly different from the one 
identified by Eq. (|4]i. These rotation curves represent systems 
with a "wrong tuning" between the disk and halo components, 
in the sense described in Sect. |2] We adopt the same simple 
parametrization for the rotation curve shape as in Eq. (l25T l and 
consider two different values for the parameter r, taken to be 
significantly far from the value which that reproduces the cor- 
relation in Eq. (|4]l. On the fast rising end, we take t/ = 0.4 
(i.e. Rn/h * 0.44), while on the opposite side of a slowly rising 
rotation curve we take t, - 2.3 (i.e. Rn/h ~ 2.53). We apply 
the self-consistent and the parametric disk-halo decomposition 
to these two anomalous rotation curves and refer to the same 
function H as in Eq. ( |29] l to quantify the quality of the fit. 

The results for the fast rising rotation curve are shown in 
Fig. IB. II and Fig. IB. 21 In the parametric decomposition method 
we find again the bimodality of the fit, as noted in the case of 
NGC 3198. The self-consistent method instead identifies a de- 
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Fig. B.2. Contours of the function H for the self-consistent de- pj ^^ ^.^^^^^^^^ ^^ ^^^ ^^^^^.^^ ^ ^^^ ^^ parametric decern- 

composition of the rotation curve dehned by Eq. i25i with .• /u j c mpriN t ^u <. <-■ j <; j u 

_^ ^ M 1 — I position (based on Eqs. (yj |2Jl) of the rotation curve dehned by 



Eq. ^ with Tj = 2.3. 



V 




p 



Fig. B.3. Disk-halo decomposition associated with the best self- 
consistent fit of the case tj - 0.4. The coding is the same as in 



composition with an important stellar disk; the fit for the rotation 
curve is illustrated in Fig. IB. 31 

The results for the slowly rising rotation curve are shown in 
Fig. lB.4l and Fig. IB. 51 The degeneracy pattern is clearly present 
in the parametric decomposition, but in this case it is less marked 
with respect to the cases studied earlier in this paper. Here, the 
parametric fit is "pushed" in the direction of an insignificant 
disk, because the observed Rn is too small and can only be as- 
cribed to a dominant dark matter halo. On the other hand, the 
self-consistent decomposition appears to be able to handle also 
this case (see Fig. IB. 61 ). 

From the point of view of the quality of the fits, it is apparent 
that both decomposition methods work best when fitting the case 
of T a; 1, that is for a rotation curve that follows the empirical 
correlation illustrated in Fig.[T] This is the aspect of conspiracy 
that was introduced and described in Sect. 12.11 In this case, the 




04 




Fig. B.5. Contours of the function H for the self-consistent de- 
composition of the rotation curve defined by Eq. ( l25T l with 

T, =2.3. 



values of the function H corresponding to the best-fitting (ff,/?) 
pairs are the smallest. In particular. Figure IB. 31 shows that the 
case of a fast rising rotation curve cannot be properly described 
without including a central concentrated mass component, such 
as a bulge; the value of the residuals S in correspondence of 
the best-fit (see Fig. lB.2l l is considerably higher if compared for 
example to the good fit of Figure [12] Similarly, for the opposite 
case of a slowly rising rotation curve, it is the parametric model 
which is quite unable to account for such a shape of the rotation 
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